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Abstract 

An exact solution of the Einstein equations for a Bianchi -I universe in the pres- 
ence of dust, stiff matter and a negative cosmological constant, generalising the 
well-known Heckmann-Schucking solution is presented. This solution describes 
a universe existing during a finite period of cosmic time, where the beginning 
and the end of its evolution are characterized by the presence of Kasner type 
cosmological singularities. 
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1. Introduction 

The construction of exact cosmological solutions of the Einstein equations 
always represents a rather attractive task. In our opinion there are now two most 
important classes of cosmological solutions. One of them is the class of exact 
solutions of Friedmann-Robertson- Walker isotropic cosmological models, which 
constitutes a basis for comparison of theoretical predictions with observations 
[1]. The anisotropic Kasner solution for the empty Bianchi-I universe Q is 
the oldest time dependent cosmological solution of the Einstein equations. Its 
importance for theoretical physics is connected with its role in the description 
of the oscillatory approach to the cosmological singularity 0], which in turn 



* Corresponding author 
Email addresses: Alexander.Kamenshchikabo.infn.it (Alexander Y. Kamenshchik), 
chiara.mingarelli@studio.unibo.it (Chiara M. F. Mingarelli) 



Preprint submitted to Physics Letters B 



September 23, 2009 



appears to be a very promising topic for study in the string and M-thcory 
context 

The Heckmann - Schucking [5] anisotropic solution for the Bianchi-I universe 
in the presence of dust-like matter constitutes a bridge between these two types 
of the cosmological solutions: in the vicinity of the cosmological singularity 
it behaves as a Kasner universe, while at the later stage of the cosmological 
evolution it behaves as an isotropic flat Friedmann Universe. The Heckmann- 
Schucking solution has the following form: for the Bianchi-I universe with the 
metric 

ds 2 = dt 2 - a 2 {t)dx 2 - b 2 {t)dy 2 - c 2 {t)dz 2 (1) 
filled with dust whose equation of state is 

P = (2) 
the functions a{i), b(t) and c(t) are given by the formulae 

b(t) = b t^(t + t ) 2 / 3 ^, 

c(t) = c t^(t + t ) 2/3 - p3 , (3) 

where the exponents p\ , pi and p^ are the well-known Kasner exponents [2, 0] 
satisfying relations 

Pl+2>2+P3 = l, (4) 

Pi+Pa+P3 = l- (5) 
Usually the Kasner exponents are arranged in such a way that 

Pi < P2 < P3- (6) 

It is easy to see that the solution ([3]) is close to the Kasner solution when 
t <C to- In the limit t S> to ah the functions a(t), b(t) and c{t) are proportional 
to i 2 / 3 , i.e. their behavior coincides with that of a flat Friedmann Universe filled 
with dust. The generalization of the Heckmann-Schucking solution for the case 
of Bianchi-I Universe filled with an arbitrary fluid with the equation of state 
p = wp, w — const, where p is the pressure and p is the energy density was 
obtained in paper 0. 

In a recent paper Q the Heckmann-Schucking solution was generalized in the 
case of the Bianchi-I Universe filled with a mixture of three perfect fluids: dust, 
stiff matter with the equation of state p = p, and a positive cosmological con- 
stant. The consideration of a positive cosmological constant in combination with 
other types of matter was motivated by the recent discovery of the phenomenon 
of cosmic acceleration Q , which makes the presence of the cosmological constant 
or some other type of matter which mimics it necessary for a realistic cosmolog- 
ical model. Another argument in favour of a positive cosmological constant is 
the fact that inflationary theories for a very early Universe [9( contain an effec- 
tive cosmological constant providing a period of a quasi-exponential expansion 
at the beginning of the cosmological evolution. 
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However, the consideration of the cosmological models with a negative cos- 
mological constant has some attractive features as well. First of all, the pres- 
ence of a negative cosmological constant is compatible with supersymmetry [1(3] ■ 
Secondly, one can formulate various cosmological models, compatible with mod- 
ern observations, where the cosmological constant has a negative sign and the 
present cosmic acceleration is driven by other types of dark energy. In such 
models, the present accelerated expansion is followed by contraction of the Uni- 



verse culminating in the Big Crunch cosmic singularity Consideration of 
the Bianchi-I Universe filled with a negative cosmological constant is therefore 
interesting because in this case, one has two singular regimes: the Big Bang 
and the Big Crunch, both characterized by strong anisotropy. In this paper we 
construct an exact cosmological solution for a Universe filled with a mixture of 
three perfect fluids: a negative cosmological constant, dust and stiff matter and 
study its properties. Section 2 is devoted to the construction of the solution, 
while the third section contains some conclusions. 



2. Bianchi-I universe in the presence of a negative cosmological con- 
stant 

We shall look for solutions of the Einstein equations for the Bianchi - I 
Universe with the metric fl} filled with three perfect fluids: dust, stiff matter 
and a negative cosmological constant, whose energy-momentum tensor has the 
form 

= diag{p, -p, -p, -p). (7) 
We represent the functions a(t), b(t) and c(i) in the following form: 

a(i) = #(*)exp(-2a(i)), 
b(t) = R(t) exp(a(i) - /3(t)), 

c(t)=R(t)exp(a(t)+0(t)), (8) 

where R(t) is the conformal factor while the functions a(t) and (3{t) characterise 
the anisotropy of the model. The components of the Ricci tensor have the 
following form: 

R = ~ |3^ + 6a 2 + 2/3 2 ) , (9) 
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Using the isotropy of the energy-momentum tensor ([7]) one has 

R\ = Rj = R 3 3 . (13) 
From this equation it is easy to obtain equations for the functions a(t) and /3(t): 

& + 3-a = 0, (14) 
R 

(3 + 3^ = 0. (15) 
From Eqs. (fT4")l . (|T5T) it follows immediately that 



(16) 



0=§, (17) 

where ao and /3o are some positive constants. 

The 00 component of the Einstein equations now has the form 

R 2 .o P 2 4ttG 

je =aa + T + — * (18) 

Choosing a convenient normalization of the constants one can represent Eq. 
(|18p for a Universe filled with our mixture of three fluids in the following form: 

R 2 -2 P 2 a M S 

R~ 2=a + T- A+ V + R^ (19) 

where —A is a negative cosmological term. For convenience we use A > 0, while 
the constants M and S characterise the quantity of dust and of stiff matter in 
the Universe, respectively. 

After a substitution into Eq. (fH?)) the expressions for a and from Eqs. 
(fl"6|) and (| 1 T[) we come to the following equation for the conformal factor R(t): 

R 2 . Af Sn /„„s 

T 2 =- A+ V + ^ (20) 

S = S + a 2 + ^. (21) 

Notice that the coefficient So includes contribution of both the anisotropy and 
of the presence of stiff matter. The fact the effective contribution of anisotropy 
to the effective Friedmann equation for the conformal factor behaves like that of 



where 



stiff matter was first discussed in Ref. [12j . The result of the explicit integration 
of Eq. (HI is 

R 3 (t) = g + ^sin(3x/Ai) - ^cos(3VAt). (22) 
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Here the integration constant is chosen in such a way to provide the initial 
condition 

i? 3 (0) = 0. (23) 
It is convenient to also have another form for the solution 



M V4S"oA + AP 



2A 



2A 



sVKt — arcsin 



M 



v / 4S* A + M 2 



(24) 



It is seen from Eq. (|2"4")1 that the conformal factor R(t), which is equal to zero 
at the initial moment (Big Bang) expands until some maximum value and then 
begins contracting, vanishing at some moment tsc (Big Crunch singularity). 
The value of tsc is 



tBC 



1 



3\/A 



7r + 2arcsin 



M 



V45 A + AP 



(25) 



Now, looking at formula (|25|) we see that the time of existence of the Uni- 
verse increases with the growth of the quantity of dust, characterized by the 
constant M and it diminishes with the growth of the quantity of the effective 
stiff matter So. However, changing both these coefficients one can have variation 
of the time tsc only up to factor 2. The dependence of tsc on the negative 
cosmological constant is much more significant. The time of existence of the 
Universe decreases with the growth of the absolute value of the cosmological 
constant and can acquire any value from to oo. 

Substituting the expression (f22|) into Eq. (fT6|) we can integrate the last 
equation to obtain an explicit expression for the anisotropy function a(t): 



a(t) 



In 



2y^ 3%/At 
3Ai Q &m 2 



(26) 



Here the integration constant is chosen in such a way to provide a Kasner-type 
behaviour of the function a(t) in the neighbourhood of the Big Bang singularity. 
Indeed, one can see that 



lim a(t) 



no 



t 



3y Sq ta 



(27) 



The constant t a is introduced into the expression (|26p to make the argument 
of the logarithm dimensionless. The expression for the anistropy function (3{t) 
has an analogous structure: 



00 
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(28) 
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Substituting the expressions 

_ 2-0 

/So 



a(t) 
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2 sin 



and (f2"2"| into Eqs. we obtain 
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(29) 



From Eqs. (j2"9"|) it follows immediately that at small values of t the scale 
factors a, b and c behave as 



a(t) ~ 

6(t)~i» 

c(t)~t*», 



(30) 



where 



1 
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(31) 



It is easy to check that these Kasner indices satisfy the relation (j4|) , while instead 
of the relation J5l) we obtain 



pI+pI+pI = i - g 2 , 



25 



35 3 (g 



25 



ft 



(32) 



where < q 2 < | . The relation ([32]) for the Kasner indices in the presence of stiff 
matter was obtained in [12j | . In the case q = it reduces to the standard formula 
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([5]). In this last case it is convenient to introduce the Lifshitz-Khalatnikov 
parameter u (l3j such that 



Pi 

Pi 



1 + u 



1 + u + u 2 
i 

r 



u(l + u) 

P3 = 1 , - | „,, , (33) 



where u > 1. 

Now, let us consider the behaviour of functions a, b and c in the vicinity of 
the Big Crunch singularity, i.e. when t — ► tsc- It follows from Eqs. |29|) that 
in this limit 



a(t)r 
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(34) 

It is easy to see that after the transition from a Kasner regime at the begin- 
ning of the cosmological evolution (t — > 0) to that at the end of the evolution, 
the axes, corresponding to the Kasner exponents p\ and P3 exchange their roles. 
At the beginning of the evolution the exponent pi is the smallest of three ex- 
ponents (see Eq. ([6]) ) and at the end of evolution the corresponding exponent 
(f ~ Pi) is the biggest one. Thus, it is convenient to rewrite the relations (|29p 
in the following form: 



a(t)r 


- (tBC 


-t) P 3 




^ (tBC 


-ty* 




^ (tBC 


- t) p 'l 



(35) 

Here we have introduced a new set of primed Kasner exponents, satisfying the 
relations (|4]) and (|32|) and they are ordered in such manner that p[ < p' 2 < p^.For 
the case when the stiff matter is absent (q = 0) we can parametrize the new 
Kasner exponents by the primed Lifshitz-Khalatnikov parameter u' such that 

v' U ' v 1 1 + U ' v 1 U ' (1 + U,) (36) 



It is easy to find the transformation law for the parameter u: 

, u + 2 



(37) 



Here, it is curious to remember that in the process of the oscillatory approach 
of the Bianchi-IX or Bianchi-VIII Universes to the singularity [H, Oj there are 
two types of transitions: the change of a Kasner epoch, i.e. the change of the 
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roles of the scale functions, corresponding to the Kasner exponents p\ and p 2 
combined with the shift vl = u — 1 and the change of a Kasner era, when the 
axes, corresponding to the exponents P2 and P3 swap roles and the parameter 
u is transformed into v! = — . 

u 

Now, we would like to study the behaviour of the functions a(t), b{t) and c(t) 
during the whole period < t < tsc- We shall dwell on the most anisotropic 
case when the stiff matter is absent. The function ait) begins its evolution 
contracting and finishes it also in the contraction phase. Thus, it can have 
either two extrema (minimum and maximum values) or none. The function c(t) 
increases both at the beginning and at the end of the cosmological evolution, 
hence, it also has either two or no extrema. The function b(t) increases in 
the vicinity of the initial singularity and decreases in the vicinity of the final 
singularity. It has one maximum value. The extremum condition, which follows 
from Eqs. (f29|) can be written down as 

sin SvAti ext + arcsm— ,, = = (1 — ipj) ■ ,. , % — 1,2,6. 

(38) 

For the function bit) the exponent p 2 satisfies the condition < P2 < \, 
the absolute value of the expression on the right-hand side of Eq. (|38l) is less 
than or equal to one, and, hence, this equation has always solutions. The only 
solution such that ti ext lies in the interval (0, tsc) is the maximum point 

*2 max = — t= 7T + arcsm (3p 2 - 1)— ====== - arcsm -= = . 

(39) 

For the function a{t) the extremal solutions exist if 



(1-3 P1 ) ; V * K3 °=<1. (40) 
V VM 2 + 4A^o " 

Remembering that — -| < p% < 0, one can show that the equation for the 
extremum has solutions at any acceptable value of p\ provided 

M 2 > 12A5 . (41) 

The last condition (|41[) is quite natural. The appearance of extrema shows 
that there is a period of isotropization in the evolution of a Universe. In turn, 
the isotropization is favored by the presence of dust (a large value of M) and 
disfavored by the presence of a negative cosmological constant (the greater is the 
value of A, the shorter the time of evolution tsc) and by the initial anisotropy 
of the Universe characterized by the value of the parameter Sq. If the relation 
(|4"Tj) is broken, then the extrema still exist for the values of the exponents pi 
satisfying the relation (|40[) . These solutions give the time values corresponding 
to the minimum and maximum values of the function a(t): 
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(1 - 3pi) 
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yJM 2 + 4AS 



(42) 
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^1 max 
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(1 - 3pi) 



x/M 2 + AAS 



— arcsm - 



x/M 2 + 4A5 



(43) 

Similarly the function c(t) has extrema at any acceptable value of the Kasner 
exponent | < p 3 < 1 provided the condition (|4*Tj) is satisfied. In the opposite 
case the corresponding equation for extrema has solutions for the values of the 
Kasner exponent P3 such that 



(3p 3 - 1) 



x/M 2 + 4AS 



< 1. 



(44) 



These solutions are 
1 



3VA 
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371 
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(45) 



x/M 2 + AAS a 



arcsm - 



V4A5^ 



x/M 2 + 4ASc 



(46) 



3. Conclusion 

We have found an exact solution for a Bianchi-I Universe, filled with a mix- 
ture of three perfect fluids: dust, stiff matter and a negative cosmological con- 
stant. This solution describes the universe existing during a finite period of 
cosmic time, beginning and ending its evolution in a cosmological singularity. 
In the vicinity of both of these singularities (Big Bang and Big Crunch) the uni- 
verse under consideration finds itself in anisotropic Kasner regimes, while during 
its intermediate evolution it undergoes a kind of isotropization process. The de- 
gree of isotropization increases with the growth of the presence of dust in the 
universe and decreases with the increasing of the absolute value of the negative 
cosmological constant. In our opinion, the most interesting feature of this new 
solution is the fact that it represents the transition from one anisotropic singu- 
larity to another. Thus, the generalized Heckmann-Schucking solution studied 
here can be interpreted as a very simplified model of a Bianchi-IX Universe, hav- 
ing chaotic oscillatory regimes at the beginning and at the end of its evolution. 
It would be very interesting to find the relations connecting characteristics of 
these two regimes, but it could prove to be a rather hard task. To our knowledge 
the process of evolution away from the singularity with the corresponding tran- 
sitions from one Kasner regime to another has only been studied numerically 
(see, e.g. 
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